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Abstract 

Information theoretic Broadcast Channels (BC) and MuUiple Access Channels (MAC) enable a single node to 
transmit data simultaneously to multiple nodes, and multiple nodes to transmit data simultaneously to a single node 
respectively. In this paper, we address the problem of link scheduling in multi-hop wireless networks containing nodes 
with BC and MAC capabilities. We first propose an interference model that extends protocol interference models, 
originally designed for point to point channels, to include the possibility of BC and MAC. Due to the high complexity 
of optimal link schedulers, we introduce the Multiuser Greedy Maximum Weight algorithm for link scheduling in 
multi-hop wireless networks containing BCs and MACs. Given a network graph, we develop new local pooling 
conditions and show that the performance of our algorithm can be fully characterized using the associated parameter, 
the multiuser local pooling factor. We provide examples of some network graphs, on which we apply local pooling 
conditions and derive the multiuser local pooling factor. We prove optimality of our algorithm in tree networks and 
show that the exploitation of BCs and MACs improve the throughput performance considerably in multi-hop wireless 
networks. 

I. Introduction 

The link scheduling problem for multi-hop wireless networks has received significant attention in the past few 
years lITl- lfTOl . The common assumption in these studies is point-to-point communication, that is, the possibilities 
of network information theory have not been incorporated. In this papeiQ we expand the scope of link schedulers to 
include multi-user communication scenarios using techniques developed in multi-user information theory. We first 
propose a generalized interference model to allow for such multi-user communication scenarios. We then introduce 
the Multiuser Greedy Maximum Weight (MGMW) scheduler for the proposed interference model and analyze its 
performance for arbitrary network graphs. For that purpose, we derive special conditions, that we shall call multiuser 
local pooling conditions. 

In a wireless network with shared spectrum, in general, interference prevents all point-to-point nodes from being 
used at full capacity at the same time. The general objective of the scheduling problem is to determine which links 
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to activate simultaneously in a network. A scheduling policy is said to be throughput optimal if it can keep all the 
queues stable under any stabilizable arrival rate vector, that is, any arrival rate vector for which the network can be 
stabilized. In a multi-hop network with multiple flows and fixed link capacities, the throughput optimal scheduling 
problem was initially proposed in lfT9l . The complexity of the optimal scheme, however, is very high, making 
it highly impractical to implement. Recently, researchers have focused on certain classes of network interference 
models which impose constraints on the set of links that can be simultaneously active in a network. One such 
model is the node exclusive interference model under which a node cannot simultaneously transmit and receive, 
and also cannot communicate simultaneously with more than one node. An optimal schedule for the node exclusive 
interference model, known as Maximum Weighted Matching, has a complexity, 0{N^) 1 16|, where N is the number 
of nodes in the network. A more general interference model is the fc— hop interference model, with k being the the 
minimum number of hops between any two active links (when fc = 1, we end up with node exclusive interference 
model). Maximum weight matching is NP-hard for k > 2 ||9]|. 

To address the complexity issue, low complexity suboptimal algorithms like greedy maximal scheduling have been 
proposed. An example of greedy scheduling is Greedy Maximal Matching (GMM) for node exclusive interference 
models H]. Apart from being suitable for distributed implementation lfT2]| . GMM has the property that at each time 
slot the sum of the weight of the scheduled links is no less than a fraction 1/2 of the maximum weight IHl, ifTSl . 
This also leads to the conclusion that it achieves at least a fraction 1/2 of the capacity region of the network [1]. 
However, the performance of the GMM scheme turns out to be far better than this lower bound in many scenarios, 
as shown in lfT4l and |l6|. The authors in |6| characterized the performance of the GMM scheme using a parameter 
called the local pooling factor, which is obtained from the knowledge of the network topology, and interference 
constraints. It was shown using this local pooUng factor that GMM was in fact throughput optimal for many classes 
of network graphs including all tree networks, under the node exclusive interference model (l3|. 

The past work on scheduling mainly focused on orthogonal resource sharing, i.e., if a link is active no other 
interfering link can be active simultaneously. Link models arising from the development of network or multi- 
user information theory have not been incorporated. For example, using superposition coding, a node could 
simultaneously transmit to two or more links at a rate lower than the individual link capacities, but higher than 
what could be achieved by time sharing between the individual links ifTSl . ifTTl . Similarly, by using successive 
interference cancellation techniques at the receiver node in a Multiple Access Channel, two or more nodes could 
transmit simultaneously to a receiver node with the achievable rate region being larger than the time sharing region. 
In a network, nodes may form information-theoretic broadcast and multiple access channels using the appropriate 
multiuser encoding technique to exploit the entire capacity region of the associated multiuser channels. In this work, 
we design a multiuser greedy scheduling algorithm, MGMW, for networks with multiuser channels. In order to 
analyze the performance of MGMW, we develop an interference model and certain associated conditions, which 
we refer to as CTjj-local pooling conditions. These new conditions involve the rates achievable over the multiuser 
channels and are different from the classical local pooling conditions developed for the point-to-point paradigm. 
Based on the multiuser local pooling conditions, we derive a multiuser local pooling factor, a'^j, and show that 
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(a) Broadcast capacity region. (b) MAC capacity region. 

Figure 1: The capacity regions of the two-user Gaussian BC and MAC are illustrated in this figure. 

the performance of our algorithm can be characterized using crf^. Indeed, we show that the efficiency ratio of 
MGMW, defined as the largest fraction of the network stability region that can be stabilized by MGMW cannot 
be less than a^,. We also show that the efficiency ratio of MGMW can not be larger than a",, another parameter 
derived from the CM-local pooling conditions. To illustrate our ideas clearly, we focus on networks with 2-user 
multiuser channels in this paper While the generalization to the n-user case is straightforward, the corresponding 
local pooling conditions are mathematically cumbersome; consequently we do not present them in this paper. 

We start the development of MGMW and its performance in Sections |III] and |IV] respectively, for multiuser 
channels in which the rate is a given fixed point on the multiuser capacity region. We then generalize MGMW 
and analyze its performance in Sections |V] and [VI] for the case in which the rate of multiuser channels is variable, 
chosen appropriately over the entire capacity region. A surprising conclusion we deduced from our results is that, 
for certain network configurations, variable-rate MGMW can lead to an inferior performance compared to fixed-rate 
MGMW. For both fixed-rate and variable rate MGMW, we consider examples with different network topologies to 
understand what topologies lead to throughput optimality for MGMW. If MGMW is not throughput optimal, we 
compare its performance with the optimal performance as well as the performance of GMM, the greedy scheduler 
without multiuser channels. We also simulate MGMW using an arbitrary network topology in Section IVlIl 

II. SYSTEM MODEL 

We model the wireless network as a graph G{V, £) where V is the vertex set representing the nodes and £ is the 
set of edges. Each edge represents a directed point-to-point link over which a sender node transmits to a receiver 
node. To present ideas clearly, we assume a single-hop network traffic model. The extension to multi-hop traffic 
models is readily possible using techniques in Q, ||4l. We assume a time slotted model indexed by t in which 
the packets arrive at the start of every time slot. The packet arrivals in each link are independent and identically 
distributed across time slots, but may be correlated across Unks. Each node keeps a separate queue for every edge 
it transmits over. Let Q, (t) represent the queue length for the packets to be transmitted over edge I and A, denote 
the arrival rate for edge I. 

We assume that a node can communicate to a single node or a pair of nodes simultaneously, using information 
theoretic broadcast and multiple access channels. The information theoretic Broadcast Channel (BC) occurs when 
a transmitting node sends jointly encoded data using a suitable codebook, intended to more than one receiver at the 
same time, from which the receivers can decode their respective information. Similarly, a Multiple- Access Channel 



4 



(MAC) occurs when a node receives simultaneously from more than one transmitter nodes . The MAC (or BC) 
capacity region for a given power constraint is the closure of the convex hull of the set of rate vectors such that 
there exists codebooks at this rate, with an average power below the power constraint, such that the receivers can 
decode with arbitrarily small probability of error ifTTl . Fig. [T] illustrates the capacity region of a two user AWGN 
Broadcast Channel and an AWGN MAC. Given the channels, and given certain coding schemes, we have certain 
BC and MAC rate regions that our rates belong to. In the rest of the development, for the ease of exposition, we 
shall refer to these as capacity regions. The important observation is that the MAC and BC rate regions are strictly 
larger than the respective time-sharing achievable rate regions. Any point in the interior of the capacity region is 
achievable by choosing the appropriate set of codes. In this paper, we restrict attention to BCs and MACs comprised 
of only two edges, i.e., a sender can transmit to two receivers (BC) or two senders can simultaneously transmit to 
one receiver (MAC). In the sequel, we use the term link in a generic sense to include point-to-point links as well as 
links made up of BC and MAC links. A point-to-point link consists of a sender node transmitting to a receiver node 
over one edge. A multiuser link is formed when a node transmits to two receiver nodes over a pair of edges (BC), 
or when a pair of nodes transmit to a common receiver node over two edges (MAC). We will specify the nature of 
the link wherever necessary. Note that the generalization of our results to the scenario in which a multiuser channel 
can involve more than two edges is straightforward. However, the algebra involved and the conditions we derive 
become highly complicated and we do not present this generalization in this paper. 

We assume that a multiuser link {k, I) can operate at any rate on the boundary or the interior of the capacity region, 
denoted by TZ^i, and that the scheduling policy has the freedom to choose the appropriate rate point (Cfc,(t), Cik{t)) 
in each time slot. In Section |III1 we initially assume a simpler network model where each transmitter, capable of 
multiuser communication chooses a single rate pair from the boundary of the capacity region for the associated 
multiuser link (BC or MAC), and whenever it chooses to transmit over that link, it uses the associated rate pair. 
Thus, for a multiuser link {k,l), {Cki{t),c,k{t)) = {Cki,Cik) is fixed throughout and the scheduler does not have 
the freedom to choose the rate of the multiuser link. In Fig. [l] (Cfc,,Cifc) is the fixed rate pair associated with the 
multiuser link (fc, I). In Section|V] we relax this assumption so that the multiuser link can operate at any point in the 
capacity region, with the scheduling policy having the freedom to choose an appropriate rate point every time slot. 
Note that the point-to-point link capacities also belong to the multiuser link capacity region (as (c^, 0) and (0, c,)), 
but are defined separately to distinguish these from the multiuser link case, where both edges of the multiuser 
link have non-zero transmission rates. For a Gaussian BC or MAC, as seen in Fig. [T] c^i < c^, and Ci^ < Ci and 
7^ + 7^ > 1. This means that we can achieve rates strictly larger than those achieved by time sharing between the 
two point-to-point links. It is worth stressing that using a BC link does not mean that the same message is being 
sent to both receivers, but rather different messages (packets) are sent to the receivers connected by the edges I 
and k. Similarly, using a multiple-access channel entails the two transmitters sending different information to the 
common receiver. It should also be noted that we do not rule out the possibility of any edge being utilized as a 
point-to-point link, even if it is a part of some multiuser link. 

To incorporate the possibility of such information theoretic BC and MAC links, we are motivated to introduce a 
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generalized binary interference model. Similar to classical binary interference models, each edge Z e 5 is associated 
with a set consisting of all edges that conflict with I, i.e., the links that absolutely cannot be scheduled when edge I 
is scheduled. We call this set the main interference set and denote it by X,. For edge /, let denote the set of edges 
that can be paired with I to form a multiuser link. We call the secondary interference set of I, to distinguish it 
from the main interference set, and also due to the fact that two edges that lie in each other's secondary interference 
set do not necessarily exclude each other: they simply reduce each other's rates. Note that and X, are mutually 
exclusive sets and if edge fc G YJ, then Z £ Y^. Let c, be the individual capacity of the point-to-point link I. Link I 
can be active at rate c, only if no other edge A: G X, U YJ is active. If edge A: G YJ is active simultaneously as edge 
I, then it implies that they are active as multiuser link (fc, /), at some rate (0^, (i), (i)), chosen from the boundary 
of the capacity region, as illustrated in Fig. [T] We also observe that the notion of main and secondary interference 
sets could serve a more general purpose than allowing for multiuser Unks. For instance, it is possible to extend the 
definition of secondary interference sets to include Interference Channels ll20l . a scenario in which two or more 
interfering links can be active simultaneously at reduced rates using suitable coding techniques. We use the term 
interfering links to indicate that these links are not allowed to be simultaneously active in a protocol interference 
model. Our interference model, while incorporating multiuser links, does carry over some limitations of the protocol 
interference model, namely the discrete nature of interference. The rate of a link ideally depends on the interference 
caused by other links, and this is captured in more realistic models like the SINR model. Scheduling problem with 
the SINR model has also been investigated extensively in the literature [10]-[15]. However, the problem of optimal 
scheduling with the SINR model is known to be NP hard and this fact is one of the main motivations of the use 
of approximate discrete interference models such as the binary interference model. Owing to a significant portion 
of the Uterature focusing on such models, there is a better understanding of wireless scheduling for graph based 
or binary interference models. We leverage this understanding in our extension of the protocol interference model 
when we introduce the idea of secondary interference sets. Further, the capacity regions of the multiuser links in 
the presence of interfering links may also not be known. However, our model only requires knowing an achievable 
rate region for the multiuser link which is strictly convex. 

We define a rate allocation vector 1*^1 of link rates where r, represents the rate of transmission over the edge I. 
A rate allocation vector must satisfy the following constraints: 

(i) If r, > then =0, \/k ^ Xi. This condition describes the main interference constraint for a point-to-point 
link I. 

(ii) If r, > and > 0, and also if fc G YJ and I G Y^, then r, = c,fc, and = c^i. Furthermore, = for all 
j G Yi U Yfc where j ^ k,l . This condition captures the constraints arising from the secondary interference set 
Yi: If a multiuser Unk is scheduled then the edges belonging to the secondary interference sets of either of the 
two edges that constitute the multiuser link cannot be scheduled. Thus a node is allowed to transmit or receive 
simultaneously over at most two edges. 

(iii) There exists no j G £ such that j does not interfere with any link and yet is not scheduled. 
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Figure 2: Five edge network with point-to-point link capacities Ci, C2, C3, C4 and Cg. Links (1,2) and (4,5) each form 
broadcast links. 

Let R denote the set of all possible rate allocation vectors r on a network graph. In general, R can be an uncountable 
set, since it includes rate allocation vectors corresponding to every rate pair on the boundary of the capacity region 
of a multiuser link. For any subset E C £, Re is defined as the set of all -dimensional rate allocation vectors 
that satisfy the constraints in (i), (ii) and (iii), with E substituted for £ in constraint (iii). As an example, a 
network graph comprised of five edges is shown in Fig. |2] along with the set of feasible rate vectors. Table H] 
describes the interference sets that we could define for this network under a node-exclusive interference model. In 
the figure, node ds can set up a multiuser link to send data to {di, d^) at a rate (ci2, C21), chosen from the capacity 
region. Similarly, node o?4 can transmit to (da, ds). The rate allocation vectors r for this network are: [0 Ca C5], 
[0 C3 C4 0], [ci C5], [ci C4 0], [ci2 C21 C5], and [ci C45 C54], where (cia, C21) and (045,054) are non zero 
rate pairs chosen from their respective multiuser capacity regions denoted by 7^12 and 7^,45. Note that in the absence 
of the BC links only the first four rate vectors would be available. 

The optimal stability region, or equivalently, the capacity region of a network is the set of all arrival rate vectors 
such that for any arrival vector in this set, there exists some scheduling scheme that can keep the queue lengths from 
growing unbounded. Here, we use the term optimal stability region to distinguish the network capacity region from 
the multiuser information theoretic capacity region. The optimal stability region of the network |19| is given by the 
interior of the set A = {A : A ^ 0, for some e Co(R)}, where Co(R) denotes the convex hull of the vectors in 
R and < represents componentwise inequality. Let tt : Q(t) — > R be a scheduling policy that selects a feasible rate 
vector for every time slot, based on the queue length state vector Q{t). Let 11 denote the set of all such scheduling 
schemes tt. For this model, the entire capacity region can be achieved by the Maximum Weight scheduler |fT9l . 
which at every time slot t, selects the rate vector which has the highest sum of queue-weighted rates. To compare 
the advantage of using multiuser links, we also define a set of scheduling policies that cannot utilize the possibility 
of multiuser links. Let R denote the set of rate vectors for a network with sole point-to-point communication. Set R 
satisfies the following interference constraints: r, > => rj. = 0, Vfc G AT, U Yi. Let tt : Q{t) — > R be a scheduling 
scheme designed based on this constraint and let II denote the set of all such schemes tt. Note that 11 C 11. Our 
objective is to find a low complexity scheme that belongs to the set 11 and characterize its performance with respect 
to the capacity region, as well as to compare its performance to that of other schemes chosen solely from 11. In 
the next section, we describe the MGMW scheme for the network model in which the multiuser link rates are a 
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Table I; Interference sets for the five edge network of Fig. ID 



fixed point on the boundary of the capacity region. 

III. Multiuser Greedy Maximum Weight (MGMW) algorithm 

Let the rate of any multiuser link {k,l) be a fixed point {Cki,Cik) on the boundary of the capacity region so 
that Cki{t) — Cki and Cik{t) — Ci^. In this case, the set of all rate allocation vectors R is now a finite set. For the 
network model with fixed multiuser link rates, we present a "greedy" scheduling policy, MGMW which selects 
a rate allocation vector from R in each time slot. MGMW, in principle is similar to the GMM, discussed in 16J. 
Before giving a precise definition of MGMW, it will be instructive to summarize its operation descriptively: 

Each link is assigned a weight, which is basically the queue-weighted link rates. In each time slot, MGMW first 
greedily picks the link (point-to-point or multiuser) with the highest weight. It then removes all interfering links 
and picks the link with the highest weight from the remaining links. This process goes on until there are no more 
links left to pick. More precisely, let C denote the set of all links (point-to-point as well as multiuser), i.e., 

C^{£yj {{k, /) G r I fc e 1^ and ? e y,}} . 

For any element m € C, we define the weight of a link W^{t) as follows: 

{QJt)ci, m is a point-to-point link j 

F F _ (1) 

Q.(i)c., + g,(t)c,,, TO is a MAC/BC link (fc, I) 
MGMW operates as follows. At any point in the algorithm, let Z denote the set of currently unselected Unks that 
do not interfere with any of the selected links. MGMW initializes Z to C and repeats steps 1 — 2 until Z = 0. 

1. Select a link m with the highest weight in Z. 

TO G argmax{PF„(t)}. (2) 

Note that to need not be unique. In case of a tie between a point-to-point link and a multiuser link, MGMW 
gives priority to the point-to-point link. 

2. After the selection, remove all links that conflict with to, i.e., set their rates in the rate allocation vector to 
zero. If TO is a point-to-point link j then the scheduler sets r(fc) = 0, for all fc G U Y^. If m is a multiuser 
link (fc, I) then it sets r{i) = 0, for all i G U U Xi L) Yi except i ~ k,l. Update Z to consist of only 
non-interfering links. 
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At the end of the procedure MGMW yields a rate vector that belongs to the set R. Also, if Yi = 0, V/ G £, MGMW 
reduces to the GMM of ||6l. 

Example 1. . MGMW for the five edge network of Fig. |2l 

The set of point-to-point, and multiuser links for this case is given by C = {1,2, 3, 4, 5, (1,2), (4, 5)}. Let the 
link rates be Ci = 4, Ca = 6, c., = 2, C4 = 8, C5 = 5, (045,054) = (4,3) and (ci2,C2i) = (3,4). Let Qx{t) = 20, 
Q^{t) = 5, Qzit) = 2, Qi{t) = 12 and Q.^i) = 1- Applying the MGMW algorithm, hnk 4 is observed to have the 
highest weight of 96. Link 1 as well as (1, 2) each have weight 80. Link 4, having the highest weight, is picked first 
and following step 2, the interfering edges 2 and 5 given in Table |T] are removed from set Z. Among the remaining 
links, highest weight is seen to be 80 for link 1. Node is hence selected to transmit over link 1. The chosen rate 
allocation vector is then [40080]. Thus, at time t, no multiuser link is chosen to transmit. 

A. Performance Characterization of MGMW Scheduler 

We adopt the definition of efficiency ratio given in |6| to describe the performance of the MGMW algorithm. 
The efficiency ratio of the MGMW scheduling algorithm 7 is defined as the largest fraction of the capacity region 
such that any arrival rate vector inside this region can be stabilized by MGMW, i.e.. 



We study the efficiency of the MGMW algorithm for any network by relating it to a parameter that we call the 
multiuser local pooling factor, which depends on the network topology and the interference sets. In the no multiuser 
link scenario, i.e, when the set = for all I, fT4| showed that the GMM scheduler is throughput optimal for 
network graphs which satisfy certain conditions. These conditions, known as local pooling conditions are based on 
the network topology and the link interference constraints. In [6|, a more general condition called u-local pooling 
was introduced to characterize the performance of GMM for arbitrary interference graphs, including those for which 
GMM was not throughput optimal, s In this section we identify new network conditions in the presence of multiuser 
links, which we call multiuser local pooling fo-jf -local poohng) conditions. We will use these conditions to define 
the multiuser local pooling factor for any network graph. Recall that C is the set of all links for a given network 
graph. To describe the (Tj(j-local pooling conditions, we focus on certain subsets of C, which we call candidate 
maximum weight (MW) subsets. A set of links Lmw C Eb is called a candidate MW subset, if there exist queue 
lengths, not all zero, such that Lmw = argmax,^^2{T4^„}. Not all subsets of £b can be candidate MW subsets. In 
fact, every candidate MW subset Lmw satisfies the following property]! For any pair of edges k and I such that 

^Consider any Lmw- Suppose there ai'e two point-to-point links k £ Lmw, I G Lmw such that I (HYi^. Then the weights of the Hnks k 
and I are equal, i.e, qkC^ = qici. Suppose the weight of the broadcast link is not greater than that of the individual links, qkCki +qicik < g^Cfc. 
Substituting for qi from q^Ck = qici gives us the condition + < 1- This contradicts our earlier assumption on the rates for the broadcast 
channel that -I- > 1, and hence the weight of the broadcast link exceeds that of the individual links. Thus, Lmw cannot be the set 
with the highest weight. 



7' :— sup {7 I the system is stable under MGMW 
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k E Y,, both individual edges k and I of the multiuser link (fc,Z) do not appear as point-to-point links in Lmw 
separately. In other words, if an edge j appears as a point-to-point link in the set Lmw, then no other edge in the 
secondary interference set of j appears as a point-to-point link in Lmw- {j S Lmw ^ I ^ Lmw for all I e YJ}. 
We also denote Slmw consisting of all the edg es in Lmw- 

Example 2. Consider the network graph of Fig. |2] Examples of candidate MW subsets for this graph are the sets 
{1, 3, (4, 5)} and {(1, 2), 2, 3, (4, 5)}. The sets {(1, 2), 1, 2, 3} and {1, 2, 3} however, are not candidate MW subsets, 
as the links 1 and 2 appear together as point-to-point links in both sets, while also comprising the edges of the 
multiuser link (1,2). For this network graph the set £lmw ~ {1: 2, 3, 4, 5}. 

We now introduce the idea of crj,f -local pooling conditions applied to candidate MW subsets. 

Definition 3. Let Lmw be any candidate MW subset. Then Lmw contains point-to-point and multiuser links as its 
elements. Let Rr,.,,,, denote the set of all rate allocation vectors for the set £l.,„,. Also, let Rr,.,,., C Rr.„., be 

^MW •> ^MW ' ^ MW — ^MW 

the set of rate allocation vectors on £lmw' ^^'^^ '-'^^ chosen by the MGMW policy, when links in Lmw have the 
maximum weight. Then Lmw satisfies <Tm- local pooUng if, for any given pair pi, v, where fl is a convex combination 
of the rate vectors in ^lmw '^"'^ ^ ^ convex combination of the rate vectors in ^lmw' ^'f^^'' of the following 
hold: 

(i) There exists a point-to-point link j G Lmw such that (7^/1^ < v^, or 

(ii) There exists a multiuser link {k, I) G Lmw such that UMil^kCki + IJ-iCik) < i^kCki + i^iCik- 

Condition (i) becomes the standard cr-local pooUng condition of 161, when defined for an arbitrary subset of edges 
in £. The (Tj^-local pooling condition is distinguished by the fact that it is stated only over candidate MW subsets. 
We introduce (ii) to generalize it to the case where multiple-edge links are possible, such as information theoretic 
broadcast or multiple access channels. We define a parameter, cr^j^, for a network as the supremum of all cTm such 
that every candidate MW subset Lmw of C satisfies (Tj^-local pooling, i.e., 

Um = sup {ctm I y Lmw G 'C, conditions (i) or (ii)are satisfied for every /i and i7} , 

where jl and i/ are convex combinations of the rate vectors in ^^lmw '^^'^ ^i-mw respectively. We call crf^ a multiuser 
local pooling factor. To show throughput optimality when local pooling conditions were satisfied, the authors in 
lfT4l argued that if a set of links alternately had the highest queue weighted rate in a small interval of time, and 
if they satisfied local pooling, then GMM served to bring down the highest weights in that interval. The proof 
used a fluid limit argument to find a Lyapunov function whose drift was then shown to be negative. A similar 
approach is followed in the proof of Lemma 1 in ||6l. When multiuser links are included, one needs to consider 
the weight of both point-to-point and multiuser links. This leads to local pooling conditions being defined over a 
fixed class of subsets, i.e., candidate MW subsets of links and not over all subsets of links. The reason for this 
will become evident in the proof of Theorem 1, where it is seen that while considering the set with links having 
maximum weight, one may exclude non-candidate MW subsets, as links in these sets cannot have the maximum 
weight simultaneously. 
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Based on the local pooling conditions, we identify a condition for a candidate MW subset ij^v, with which 
aM-local pooling does not hold. L^w does not satisfy Gm-IocqI pooling if there exists a pair of vectors jl, v that 
are convex combinations of the rate vectors in and Rt^^j,, respectively, such that they satisfy o-j/yU >r i'- We 

define the parameter as follows: 

a" J = mi{aM \ 3 an L^iv G -C, such that CTj//! ^ i7 for some jl, 

where z7 is a convex combination of rate vectors in Rr.„„ and a is a convex combination of rate vectors in Rr„„,. 
We now give Theorem [T| and Theorem |2] to prove that the efficiency ratio of the MGMW scheduler satisfies the 
relation cTj^ < 7* < cr",. 

Theorem 1. The network is stable under MGMW algorithm for all arrival rate vectors A satisfying A G cTj^A. 

Proof: The Proof is given in Appendix A. ■ 
While Theorem [U shows that any arrival rate within cr^jA is stabilizable by the MGMW algorithm, we further 
link the performance of the MGMW to a", in Theorem |2] by showing that there exist arrival rates, arbitrarily close 
but strictly outside of for which the system is unstable under the MGMW scheme. Theorem [1] together with 
Theorem |2] implies that the efficiency ratio of MGMW is bounded below by trfj and bounded above by (7^^. 

Theorem 2. Let there exist a candidate MW subset L^,n- G C such that for some positive number Um, ond a pair 
of vectors which are convex combinations of the elements in ^L^.ny <^nd ^L^.^y, Cm/I >n v is satisfied. Then, 
for any e > 0, there exists a k )~ such that the arrival rate \ = v + ek makes the system unstable under the 
MGMW scheme. 

Proof: The proof is given in Appendix B. ■ 
Theorem 1 and Theorem 2 establish sufficient and necessary conditions respectively for local pooling to occur under 
the MGMW policy. Unlike the GMM policy, for which a single condition is both necessary and sufficient for local 
pooling, we have separate necessary and sufficient conditions due to the presence of multiuser links. While these 
two conditions are identical in certain cases such as tree networks satisfying specific rate constraints, they are not 
necessariry identical in a general network. 

IV. Performance of MGMW 

In this section, we analyze the performance of MGMW in some sample network topologies. We use the bounds on 
efficiency ratio of MGMW obtained in the previous section to evaluate the throughput gain by leveraging multiuser 
links in these network graphs. Even though the optimal stability region of a network with multiuser links is larger 
than that of a network with point-to-point links alone, for some networks it may be possible that MGMW achieves 
a smaller stability region than GMM. However, we provide examples of two network graphs, specifically a tree 
network graph and a star network graph, for which we show that MGMW achieves a larger stability region than 
GMM. In this process, we also explore the tightness of our bounds. In particular, we show that the lower bound is 
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tight for tree networks when the multiuser link rates are restricted to certain rate regions. Moreover, the lower bound 
on efficiency ratio is good enough to prove that the stability region with MGMW is larger than that with GMM 
for the star network example. We also simulate the performance of MGMW and GMM on a randomly generated 
network graph, and observe that MGMW can stabilize a higher arrival rate for each link in the random network 
graph. 

We first show that MGMW is throughput optimal for certain tree networks. 

A. Tree Networks 

In the following theorem, we show that under the node exclusive interference assumption, MGMW is throughput 
optimal for directed tree graphs in which no two multiuser links have a link in common. The node exclusive 
interference assumption for our model only restricts a node from transmitting and receiving at the same time. 
It does not restrict a node from transmitting simultaneously to multiple nodes, or receiving simultaneously from 
multiple nodes. 

Theorem 3. Let Q — (V, £) be a directed tree graph such that Yi- H =0 for all k,l E £. Let cl^^ + cf^^ > 
max(c„,c„„, c„c„„) for every multiuser link (m, n) in C. Then, if the primary and secondary interference sets are 
constructed under the node exclusive interference assumption, cr^, = — 1, implying that MGMW is throughput 
optimal for this tree network graph. 

The proof is given in Appendix ICl Theorem[3]also shows that MGMW is throughput optimal for downlink cellular 
networks without intercell interference with multiuser links consisting of broadcast channels. This is because the 
downlink cellular model is an instance of a tree network. 

B. A Network with a", < 1. 

Consider the network graph shown in Fig. [5] All links have a rate of 1 when used as point-to-point links. All 
multiuser links have a rate of (0.75,0.75) each. We define cr^jj^,,- for a candidate MW set Lmw as the highest value 
of (TAj-local pooling satisfied by Lmw- We show in 121] that cfl^w ^ 2/3 for all candidate MW subsets Lmw of 
the star network graph in Fig. [3] and hence the efficiency ratio of MGMW is at least 2/3 for this network graph. 
Here we will only show this for one candidate MW subset and will not repeat the same operation for all candidate 
MW subsets. In order to find (Tlmw^ make use of the following relation that we derive in Appendix ICl 

^--max,,...,,_,||r^-T||,' 
where f' e ^lmw^ ^ ^i-^w T is a \£lmw \ ^ I^^-mw] matrix such that T„ = 1, if i e Lmw', T,^ = 
Cki and r„ = c,fc if (fc, /) G Lmw\ and = otherwise. 

Consider the candidate MW subset Lm„ = {(1, 2), (3, 4), (5, 6), (7, 8), (9, 10), (11, 12)}. The set R^^^^ is given 
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Figure 3: A twelve edge network with six broadcast links. 
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Applying the relation (|32] | to the specified Ljvf,,i' set yields ctj^^^^^^, ^2/3. Furthermore, we also show in [1211 that 
< 3/4 for this network graph. For the same network graph, if we exclude the possibility of multiuser links, 
it can be shown that GMM algorithm has an efficiency ratio of 2/3 (with respect to the capacity region in the no 
multiuser link scenario). To see this, we only need to consider the set L = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12}. The 
set L forms a cycle, and in the same manner as that shown for the six-link cycle network in [6|, one obtains the 
local pooling factor as 2/3. Hence, in this example, MGMW guarantees a larger stability region compared to GMM. 
This is because the optimal stability region of a network for the multi-user case is a superset of the optimal stability 
region of the same network in the no multiuser link scenario. 

In the following section, we explore a more general network model, where we assume that rather than just one 
point, the whole rate region is available to the encoder and decoder, as well as the scheduler. The scheduler then 
selects a suitable multiuser link rate pair from this region in every time slot. For this network model, we provide a 
generalized version of the MGMW scheme and analyze its performance by deriving local pooling conditions that 
are similar to the ones derived for the case with fixed multiuser link rates. 

V. Using the entire multiuser capacity region: variable rate MGMW 

We present the extension of the MGMW scheme wherein we allow it to select an arbitrary rate pair from the rate 
region of each multiuser link in each time slot. From this point on, we restrict attention to network graphs that only 
have BC links as multiuser links, and omit the technical treatment of MAC links due to space constraints. While 
the variable rate MGMW algorithm remains unchanged in the presence of MAC links, the performance analysis of 
variable rate MGMW involves treating the MAC links and BC links as separate cases. 

The MGMW scheme, as described before in Section |III1 initially determines the weight of all point-to-point and 
multiuser links. Since the entire multiuser capacity region is now available to the scheduler, the weight of a multiuser 
link is defined as the sum of queue weighted edge rates, maximized over all possible choices of rate pairs within 
the capacity region of the multiuser link. For any element m G £, We define the weight of a link W^{t) as follows: 
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Qj{t)Cj, m is a point-to-point link j 



max Qfe(i)cfci + (5i(0cifci (5) 



m is a multiuser link (fc, V), 

where i?fc, is the capacity region of the multiuser Hnk (fc, /), and is assumed to be known. For instance, consider 
the operation of variable rate MGMW scheme for the network graph of Fig. |2] 7^i2 and 7?.45 are the multiuser rate 
regions of links (1, 2) and (4, 5) and are assumed to be known. The variable rate MGMW then uses (|5j to calculate 
the weight of each link. Let (c\^(t)^ c^^it)) and (cl^it), c*j(<)) be the transmission rate pairs that correspond to the 
weight of the multiuser links (1, 2) and (4, 5). Now that all the weights and their associated transmission rates are 
known, variable rate MGMW then assigns the rate allocation vector following the same procedure as described for 
the MGMW scheme in Section |III] 

For a BC link, since the capacity region of the multiuser link is strictly convex, as shown in Fig. [T] the transmission 
rate pair associated with the weight of the link is unique, and corresponds either to a multiuser link configuration (if 
the rate pair that yields the maximum weight allocates non-zero rates to both edges) or to one of the point-to-point 
links (if the rate of one of the edges is zero). Thus, only one of the configurations can be associated with the weight 
of the multiuser link at any given time. With the new weights calculated as described above, the operation of the 
MGMW algorithm is identical to the one described in Section |III] and yields a rate allocation vector that belongs to 
the set R. The set R for a network graph is now an uncountable set, as it includes every rate pair on the boundary 
of the capacity region, rather than one fixed rate point. 

Local pooling conditions for variable rate MGMW 

The performance of the variable rate MGMW scheme can be analyzed in a manner similar to that of MGMW, by 
relating the multiuser crj,^ -local pooling factor to the efficiency ratio of variable rate MGMW. We retain the same 
notion of efficiency introduced earlier in Section |IV] Note that the network capacity region defined in Section |II] is 
a superset of the capacity region obtained by fixing the multiuser link rates. We now present the multiuser cr ji^-local 
pooling conditions for the variable rate MGMW scheme. 

Let Lmw be any candidate MW subset, i.e., a subset of C whose links can achieve the maximum weight simul- 
taneously. Every candidate MW subset L^w satisfies the following property: For any pair of edges k and I such 
that k and vice versa, a candidate MW subset L^w can contain at most one out of the three links (fc, /), k 
and /, i.e., 

for any multiuser link (fc, /), where '^n^L^^^^ indicator variable taking the value 1 when ii G Lmw 

and 

otherwise. This property of L^w follows from the fact that the weight of link (fc, I), as a consequence of being 
maximized over the entire multiuser capacity region TZ^i, can correspond to only one out of the three possible 
configurations (point-to-point link k, point-to-point link I or multiuser link (fc, I)). 

Given a candidate MW subset, the transmission rates of the point-to-point links are fixed at their respective capacity 
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values, whereas the transmission rates of multiuser links can take different values from the associated capacity 
regions. Let Ci^^j,, be the set containing all possible transmission rates associated with the links in Lmw We 
illustrate these sets in the example below: 

Example 4. Consider the network graph of Fig. |2] An example of a candidate MW subset for this graph is the 
set {1,3,(4,5)}. The sets {(1,2), 1,3} and {(1,2), 2, 3} however, are not candidate MW subsets, as they both 
contain links (links (1,2) and 1 in one, and (1,2) and 2 in the other) that cannot have the maximum weight at 
the same time. For the example Lmh- = {li3,(4, 5)}, £lh.,„, is {1,3,4,5} and the edges have rates given by 
{ci,C3, (c*5,c*^)}, where (0*5,0*4) is an arbitrary non-zero rate pair selected from the boundary of TZi^. Since the 
rates of point-to-point links are fixed, each element of Ci^^^^ would correspond to a non-zero rate pair chosen from 

7^45. 

We now state the CTjvj-local pooling conditions applied to candidate MW subsets. 
Definition 5. Let R, be the set of all rate allocation vectors on the subset of edses E, For each transmission 
rate vector c G Cr let R? be the set of all rate allocation vectors that MGMW can assisn on £, when 
the links in Lmw achieve maximum weight with their transmission rates given by c. Note that ^^l^w ^^mw ^ 
^i-MW then say that Lmw satisfies CTij-local pooling, if every ^'lmw ^'^^^^fi^^ the following conditions: 
For any two vectors fl and v, where /Z is a convex combination of elements in Rlmw' '^^'^ ^ ^ convex combination 
of rate vectors in R? 

•' ^ MW 

(i) There exists a point-to-point link j £ Lmw such that ffMl^-j < i^j, or 

(ii) There exists a multiuser link (fc, /) G Lmw such that aMfJ-k < t^k ond GmI^i < i^i- 

The multiuser local pooling factor o'm of a network is the supremum of all Om such that every candidate MW 
subset Lmw of £ satisfies crj,f-local pooling, i.e., 

(t'm = sup {(Tjvj I V Lmw G £, V c e Conditions (i) or(ii) are satisfied for every fl and i7} , (6) 

where jl and z7 are again convex combinations of the rate vectors in Rlmw '^^'^ ^^mw respectively. For the 
variable rate MGMW scheme, we state the following result that extends the performance analysis results for the 
MGMW scheme with fixed multiuser link rates. 

Theorem 4. For a graph Q with the multiuser local pooling factor (t*^ as defined in (|6]l, variable rate MGMW 
scheme stabilizes all arrival rate vectors A G CTmA, i.e., 

7* > (7) 

Proof: The proof is given in Appendix D. ■ 
We observe that for the variable rate MGMW, the multiuser local pooling conditions involve verifying the 
conditions over all possible sets of transmission rates associated with the links in Lmw This is in contrast to the 
case where the multiuser link rate is fixed, since there is only transmission rate vector associated with the links in 
an Lmw subset. We now derive an upper bound for the efficiency of variable rate MGMW can further be analyzed 
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Figure 4: A network with two multiuser links. 

as follows. Let o-j,, be defined as: 

(7 J,/ = inf {(Tj,, I 3Lmw G -C,and an associated transmission rate vector cG Ci^^^^, (8) 

such that cTj,//! > for some /2, i?} , (9) 
jl and i7 being convex combinations of the rate vectors in f^Lj^.j^ and ^1^^^ respectively. 

Theorem 5. For a graph Q with Om as defined in (|9]l, there exist arrival rate vectors A arbitrarily close to, but 
outside (Taj A, such that the network graph is unstable under variable rate MGMW. 

7* < ^M. (10) 

Proof: The proof is similar to that of the fixed multiuser link rate case and is given in Appendix E. ■ 
Theorem m and Theorem |5] together imply that ct^j < 7' < (Tm- Since the ctj,/ -local pooling conditions need 
to be satisfied by all rate assignments associated with the links that belong to an L^w subset, some of the rate 
assignments could limit ctm more than other rate assignments, leading to a lower value of efficiency for variable 
rate MGMW. 

The stability region of a network with multiuser link becomes even larger when we allow the entire achievable rate 
region for each multiuser link to be used by the scheduler However, the proposed variable rate MGMW policy, 
which is a natural extension of the fixed rate MGMW policy can perform poorly compared to the fixed rate MGMW 
policy. Indeed, the efficiency ratio of the variable rate MGMW policy is in fact bounded above by the efficiency 
ratio of a fixed rate MGMW policy (with fixed rates being chosen from the rate region of multiuser links). In the 
following section, we provide examples of network graphs in which we illustrate how the availability of the entire 
multiuser capacity rate region can bring down the efficiency of variable rate MGMW. We observe further that GMM 
is actually throughput optimal for these network graphs examples in the absence of multiuser links. These examples 
illustrate that it might be more beneficial to use a fixed rate MGMW scheduler whose rates have been carefully 
chosen from the corresponding rate regions of multiuser links compared to the variable rate MGMW scheduler. 

VI. Performance of variable rate MGMW 
In this section, we illustrate the performance of the variable rate MGMW scheme using some examples. Specif- 
ically, we provide some cases in which the variable rate MGMW performs much more poorly in terms of the 
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efficiency ratio, compared to the fixed rate MGMW scheduler to the extent that, certain arrival rates that are 
stabiUzable by fixed rate MGMW make the queues blow up for variable rate MGMW. While the optimal network 
stability region is larger when the rates are allowed to vary, we show that the variable rate MGMW may perform 
worse relative to fixed rate MGMW with regard to efficiency ratio. 

Example 6. Consider the network graph shown in Fig. |7a| The network has four edges and is capable of using 
the multiuser links (^1,^4) and {l^jh), whose multiuser capacity regions are given by TZ^^ and TZ^^ respectively. 
We assume that all the point-to-point links have unit capacities. Let (1 — £1,64) G TZn and (ca, 1 — Ci) G 7^23 be 
two arbitrary rate pairs. We now show that the efficiency of variable rate MGMW is limited by the choice of the 
multiuser link rate pairs according to the relation: 

7<1±^. (H) 
where 5 = (€4 — ej + (ca — £1) is a positive quantity that depends on the value of the link rate pairs. 

Proof: Consider the candidate MW subset Lmw ~ {{h, L), {I2, L)} with the link rates given by (1 — Ci, 64) G 
7^14 and (ca, 1 — ei) G TZ23, where ei, 64 G (0, 1). Let z7 be a convex combination of the rate allocation vectors on 
Lmw, with the Unk rates that yield the maximum weight assumed to be (1 — £1, £4) and (ca, 1 — Ci). We choose 
to be: 

£7 = 0.5 * [1 - £1 £4] + 0.5 * [0 £2 1 " £1 0] . 
Let fl be a convex combination of the rate allocation vectors in f^L^w- choose fl as: 

/i = 0.5A' ((1 - £1) [1 1 0] + £2 [0 1 0] + £4 [0 1]) . 

Then, fl = Kv, which implies from Theorem |5] that the efficiency of variable rate MGMW cannot be greater than 
j^. For /i to be a convex combination of the rate allocation vectors, it must satisfy {).bK[{l — £1) + £3 + £4] = 1, 
implying that 

^ = (l-£i+£. + £4)/2. (12) 

Since the multiuser capacity region is strictly convex, 1 — £1 + £4 > 1, which implies that £4 > £1. Similarly, 
1 — £1 + £2 > 1, which implies that £2 > £4. Let 5 — {e^ — £1) + (£2 — £1), so that 5 > 0. Substituting the numerator 
in (fT2b with 5, we obtain, 

^ = {l + e,+5)/2. (13) 

■ 

Eq. ( fTTT i suggests that the efficiency of variable rate MGMW scheduler for the network graph in Fig. [Ta] is 
Umited by those multiuser rate points which have smaller values of £1 and 5, as they yield lower values of a,,,. In 
this example, a smaller £1 and 6 corresponds to those multiuser rate points which are close to one of the edges' 
point-to-point capacities, as shown in Fig. ITbjandlTcl 

Example 7. Consider the network graph in Fig. |5] with four edges having multiuser links (1,4) and (2,3) , where 
7?,i4 and TZ^^ are their respective multiuser capacity regions. Let (1 — £1, £4) G 7?.i4 and (1 — £2, £3) G 7^23 be two 
arbitrary rate pairs. Then, we show that the efficiency of variable rate MGMW is upper bounded as : 
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Figure 5: Four edge network with two multiuser links. 

2 - (ei + £2) 

7 <= ^ ^, where e, + e, < 1. (14) 

Proof: We choose Lj^m- = {{li, h), {I2, L)} with links having rates (1 — ei, e^) G TZi4 and (1 — 62, £3) G 7^23- 
Let z7 — 0.5 * [1 — £1 £4] + 0.5 * [0 1 — £2 £3 0] be a convex combination of rate allocation vectors in L^iv with 
the assumed link rates, and fl = 0.5 * K * {1 — e^) * [1 1 0] + 0.5 * X * (1 — £2) * [0 1 1] be a convex combination 
of rate allocation vectors in Rr,.,,,,. Using the fact that /2 is a convex combination, we obtain the relation: 

-i-(2-(£,+£,))/2. (15) 
In this case, the condition fl = Kv is satisfied if £1 + £3 = 1, and £2 + £4 = 1- The strict convexity of the multiuser 
capacity region gives us the conditions £4 > £1, and £3 > £2. Both these conditions together imply that £1 +£2 < 1- 

■ 

From Eq. ( fT4l l. we observe that increasing the quantity £1 + £2 decreases the upper bound on the efficiency of 
variable rate MGMW to close to 0.5, and hence rate pairs in which at least one of the edges get high rates, dominate 
the performance of variable rate MGMW. Note that in the absence of multiuser links, QMS is throughput optimal 
for both network graphs O. 

The above examples highlight certain scenarios where the variable rate MGMW scheduler can perform poorly in 
terms of efficiency ratio, when compared to a MGMW policy whose links operate a fixed rate chosen carefully 
from the multiuser capacity region. Choosing a fixed rate point also reduces the coding complexity of the multiuser 
links by requiring fewer number of codebooks. 

VII. Simulation of the Performance of MGMW 

In this section, we simulate the performance of MGMW with fixed multiuser link rates and GMM in a randomly 
connected network graph. Figure |6a| shows an arbitrary network graph having point-to-point as well as multiuser 
links. The multiuser (BC) links in this graph are links (1,2), (4,7), (3,8), (5,6), and (9,14). In this example we chose 
the transmission rates of the point-to-point links at random, uniformly between 3 and 10 units. The rates of the 
multiuser links are chosen to ensure the convexity of the rate region is assured. The arrival process for each edge 
is Bernoulli and we denote the arrival rate with A,. 

In Fig. |6bl we plot the total queue size (sum of queue lengths at each edge) as we increase the arrival rate in 
edges 9 and 14, from 2 to 2.5 as we keep A, = 1 for all other links. Here, the transmission rates of links 9, 14 and 
(9,14) are 6, 4 and (4,3) respectively. The graph suggests that MGMW yields a constant gain in arrival rate for the 
multi-user links as each edge of the multiuser link (9,14) is seen to sustain around 5% more traffic, as compared 
to the GMM case. 
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Figure 6; An Arbitrary network graph with 15 edges and 5 broadcast Unks. 

In Fig. |6c] we simultaneously increase the arrival rate across edges 1, 2, 9, 14, 3 and 8 while keeping the arrival 

rates at other edges fixed at 1. Here, the capacities of these links are 10, 8, 6, 4, 12 and 8, and the transmission 

rates of the BCs (1, 2), (9, 14) and (3, 8) are (9,5), (4,3) and (10,6) respectively. Similar to the previous scenario, 

the plots again show that the total queue size with MGMW is lower than that with GMM. Here, each edge of the 

multiuser links is able to sustain 10% more traffic than in the no multiuser link case. Thus, for the network in Fig. 

l6al MGMW appears to stabilize a larger range of arrival rates. 

Vlll. Conclusions 

In this paper we explored the problem of link scheduling in a setting that allows for the use of techniques from 
multi-user information theory. To this end, we proposed a modified version of the binary interference model by 
introducing the notion of a secondary interference set for each link of the network. The interference model proposed 
in this paper could be thought of in a loose sense as a hybrid of the binary interference model and the SINR model. 
Since the optimal algorithm is known to have high complexity (NP hard in many cases), we provided a suboptimal 
greedy algorithm called MGMW for our interference model. We characterized the performance of MGMW by 
deriving local pooling conditions and relating the multiuser local pooling factor to the efficiency of MGMW. For a 
network with capacity region A and a multiuser local pooling factor (Tm, we showed that MGMW stabilizes every 
arrival rate vector in A and that there exists a non-stabilizable arrival rate vector, arbitrarily close to, but strictly 
outside of ctj/A. We gave examples of certain network graphs where MGMW was throughput optimal and a graph 
where the multiuser local pooling factor is less than one. We also considered a network model where the scheduler 
has the freedom to select the multiuser link rate every time slot. We analyzed the performance of the variable rate 
MGMW scheduler and showed that the availability of the entire rate region could hurt the performance of variable 
rate MGMW. Finally, we also observed the performance of MGMW in an arbitrary graph and compared it to that 
of GMM. 
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Appendix A 
Proof of Theorem[T] 

The proof shows stability of the network under MGMW by finding a Lyapunov function and showing that it 
has a negative drift for the fluid hmit of the system. The idea of the proof is similar to the stability proof in 
||6|, which is for the scenario of no multiuser links. We assume that the arrival process for each link satisfies 
conditions for the fluid limit to exist, which is that the Strong Law of Large Numbers (SLLN) should hold for 
the arrival process. For example, SLLN holds when the packet arrivals in each queue are IID with bounded 
second moments. Here, we assume a modified arrival process, where we relax the IID assumption in the first 
time slot alone, and allow a deterministic but finite number of packets to arrive in the first time slot. Note that 
this does not affect the applicabihty of SLLN and Z]r=i ("-)/"- ^ = -^iW.p 1. Let A{t) denote the 
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cumulative arrival process and S{t) denote the cumulative service process until time slot t. For the arrival and 
service processes, we use A, (t) — A,([tJ), and Si{t) = Si{[t\). For the queue process Qi{t), we employ linear 
interpolation. We now consider a sequence of scaled queuing systems (Q"( ), S'"(-)). where we apply the 

scaling Q"{nt)/n, Ai{nt)/n, and Si{nt)/n, Vl E £ with the queue process satisfying l^,££(3r(0) — Then, 
using the techniques to establish fluid limit in [11], one can show that a fluid limit exists almost surely, i.e, for 
almost all sample paths and for any positive n oo, there exists a sub-sequence rij with Uj cx) such that 
following convergence holds uniformly over compact sets. For all I £ £, — A"^ ("j^) ^ ^it^ —S"^ (t) s,(i), and 
— Q"^ (Uji) — > qi{t). qi{t) and s,(i) are the fluid limits for the queue length processes and the service rate processes 
respectively. The fluid limit is absolutely continuous and hence the derivative of g, [t) exists almost everywhere ifTTl 
satisfying: 

qi{t) = <i (16) 
otherwise. 



dt 



where tt, (<) ~ -^{s,{t)). We now show that the largest queue weighted rate (taken over point-to-point link or 
broadcast links) of the fluid limit model always decreases under the MGMW algorithm. This allows us to define 
the Lyapunov function for the system as the maximum weight over all links and establish its drift to be negative. 
Consider the times t when the derivative -^qi{t) exists for all I e £. Let Lo{t) denote the set of links with the 
largest weight, i.e., 

Lo{t) = {™ G I w,^ = max„££ w,„}. 
Define the derivative of the weights of links in L as follows: 

-^qj{t)Cj m is a point-to-point link j 

" \ fMiy"' + m is a 

multiuser link (fc, /) 

Let L{t) denote the set of links from Lo{t), which have the maximum derivative of the weights, 

L{t) = {m e C I w^{t) — max w^{t)}. 

Then, one can find a small S such that in the interval {t, t + 6), links in L{t) will have the highest weights in that 
time interval, i.e., min,„eL(t) Wm(T) > max„££\f,(t) i()(r) for all t in {t,t + S). MGMW will select links from the 
set L{t) first in the interval {t,t + 6), since it picks the links in decreasing order of weights. If we focus on the 
links in L{t), then any rate allocation vector selected by MGMW in {t,t + S), projected on the set L{t) would 
yield a rate allocation vector that is an element of Ri,(t). This in turn implies that 7r(t), the service rate vector under 
MGMW projected on L{t) is a convex combination of the elements of Ri(f). A formal argument to show this is 
mostly identical to that in 1|6J and is omitted here. Because of the convexity condition ^ + ^ > L{t) is a 
candidate MW subset. Consider a A lying strictly within cr^^A. Since cTj^ is the local pooling factor, and L{t) is 
a candidate MW set, it follows from the multiuser local pooling conditions that there exists a link m G L{t) that 
satisfies 
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Xj < TTj{t), if m is point-to-point link j, (17) 

AfcCfc, + A,c,fc < 7rfc(i)Cfc,+7r,(t)c,fc, if m is multiuser link (fc,0- 

For any candidate MW set Lmw, let Sij,^,, be any convex combination of the elements (rate allocation vectors) of 
the set Rlmw- Since every Lj^,„ C C satisfies crfj local pooling, and A lies strictly within crf,A, the quantity 



€g = max 



is strictly positive for every Lmw- We then define the infimum of all such positive quantities over all such subsets 
Lmw and all vectors St„^^, as 

e* = inf . (18) 

and we observe that e* > 0. Hence, from the relation in ( fTTI i. there exists m G L(i) such that the following holds: 

Xj — TTjit) < — £*, if m is point-to-point link j, (19) 

(Afc - 7rfc(t))cfc, + (A, - 7r,(t))c,fc < -e*, if m is a link {k,l). 

From ( fT6] l, and the fact that all links in L{t) have the same derivative, ( fT9] l implies that < —e', Vm £ L{t). 
Hence, we observe that there exists no link in L{t) with > 0, and > 0, where — ma,x{q^ q,), if m is a 
multiuser link (fc, I). Now, we can consider the following Lyapunov function V{t) max^^c w^. For V{t) > 0, 
we have that 

— < maxw„ < -e*. (20) 

where -^V{t) — lim^^^o is the right hand derivative of V(t). This implies that the largest weight must 

decrease in the time interval {t,t + 6). Since the above inequality holds almost everywhere in t, the negative drift 
of the Lyapunov function implies that the fluid limit model of the system is stable and hence by the result in ifTTIl . 
the original system is also stable. 

Appendix B 
Proof of Theorem[2] 

We construct an arrival traffic using the approach in ||6l and show that under this traffic pattern the network is 
unstable under MGMW. Let J denote the number of possible rate allocations on the set L^w- The vector iJ, being 
a convex combination of the elements in Rr,,,,,, can be written as 

j-i 

i = 

where > for all < i < J— 1 and '^''^g lxj^ — 1. Let be a rational number which satisfies Ylllo I'^i ~ ''^i I ^ 7 
for any (5 > 0. Such a clearly exists for every w^. To enable the construction of the traffic pattern we now define 
a new vector v = Yllla ^i^- Thus one can make the vector v arbitrarily close to v. We now specify our arrival 
traffic with load v + ek, such that the system is unstable under MGMW. The arrival traffic for every queue consists 
of IID packet arrivals in each time slot, except for the first time slot alone. Packets arrive at the beginning of time 
slots. Without loss of generahty we assume that for the candidate MW subset Lmw, if {k, I) is a multiuser hnk 
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in Lmw, then k denotes the edge of the link {k,l) that can be present as a point-to-point link in Lmw- Let the 
initial queue state vector be Q„ >z. 0. Let j denote any point-to-point link belonging to L^w, and (fc, /) denote 
any multiuser link in Lmw Then the packet arrivals in the first time slot are such that Qa satisfies the following 
constraints: 

QjCj = QkCki + QiC,k = QkCk + Cfc - CfcCfc, > QiCi + c,Cik, (21) 

for all links j, {k,l) e Lmw- We show in Lemma [T] that a vector Qo h satisfying i2T[ exists. We now describe 
the statistics of our arrival traffic. Let vector r* be picked with probabihty ly^. Conditioned on being picked, one 
of two events may occur: 
1. With probability 1 — e, 

(a) ri{j) packets arrive into every point-to-point link j e Lmw such that ri{j) = Cj > 0, and 

(b) Cfc, and Cj^ packets arrive into links k and / respectively V(fc,Z) G Lmw such that ri{k) > and ri{l) > 0. 
We show that when packets arrive in the manner described in event 1, MGMW picks rate allocation vector r^, and 
at the end of time slot t, the queues in Lmw continue to satisfy the relation in i2T[ . 

(i) The weight of any point-to-point link j G Lmw satisying r^Q) — c^, is given by QjC^ + Using ( 1211 1 we 
obtain: 

QjCj + > QiQ, V/ G Lmw such that Tiil) = 0, and 
Q,c^ + d > Q^c„^„ + (5„c„,„ > (9„c„, 
V(m, n) G Lmw such that ri{ra) = ri{n) — 0. (22) 

The inequality in (l22T l follows since cj^ — c„,c„„, > 0. 

(ii) For any multiuser link {k,l) satisfying ri{k) — c^i and ri{k) = c^i, 

{Qk + C^i)Cu + {Qi + C,fc)c,fc > QjC,, Vj G Lmw s. t r,{j) = 0, 
QkCki + QiCik + + c^^ > (5„c,„„ + (5„c„,„, 

V(m, n) G ijv/vr such that ri{ra) — ri{n) = 0, 
QkC-ki + QiCik + cji + c'^ > QiC; + CiCik, 

QkCki + QiCik + cl, + c^^, > QfcCfc + CfcCfc,, (23) 
where (l23T l holds since 

Qk^ki Qi^ik + = Qk^-k ^~ ^k^ki (c^ c^/) -i- c^^. 

(iii) For link k G ijfvr satisfying ri(fc) = c^. and (fc, /) G Lmw for some /, we have: 

QkCk +cl> Q^Cj, G Lmw with r,(j) = 
QfcCfc + 4 > Q„c„,„ + Q„c„„ > Q^c^, 
QkCk +cl> QiCi, and 

QfcC, +c', = (Q, +c,)c,, +g,Q,. (24) 
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The weight of link k therefore dominates all the links not served by f*. MGMW breaks the tie between the 
weights of k and (fc, I) shown in (|24] | by scheduling link fc. (i), (ii) and (iii) establish that weight of links served 
in r* dominate the weight of other links and hence MGMW schedules r^. 
2. With probability e, packets arrive into the queues in L,,hy in the following manner. 

(a) Cj + Cj packets arrive into point-to-point links j G L^n- for which ri{j) — c^. 

(b) Cfc, + Cfc packets into link fc and c,fc + c, packets arrive into the queues of fc and I for multiuser links (fc, I) such 
that r;(fc) > and r;(/) > 0. 

(c) c„ packets arrive into all other links m G ij,,,,- for which r;(m) — 0. 

The quantities c^, c^, and c, are such that they satisfy the following weight criteria: 

CjCj — Cj^Ck — c^Cki CiCif. > c^C;. (25) 
It can be shown using an argument identical to that used for ( 1211 1 that there exist positive c^, Cj., and c, that satisfy 
relation (IZST l. When packets arrive according to the event 2, MGMW still schedules ff. This is because (ISTT i and 
(Us) yield 

(Q, + Cj)cj. = (Qfc + Cfc)cfc, + (Q, + c,)c,fc 

= {Qk + Cfc)Cfc + - CfcCfci 
> (Qi + C,)Ci + CiCi^. 

which again satisfies relation (1211 1. However, at the end of time slot t, when packets arrive as in event 2, the length 
of each queue j G ^tjj^„ increases by a fixed positive quantity Cj. We can now describe the queue evolution for 
our arrival traffic. The initial queue state satisfies dSTT i. Therefore at the end of each time slot, with probability 
1 — e, the queues of all edges j e £l,^.,„ remain unchanged and with probability e, the queues increase by a fixed 
positive quantity. Since the queues in Lj^^m, are non-decreasing, and the event that the queue length increases by a 
fixed positive quantity occurs infinitely often, the system is unstable under the MGMW scheme. The arrival rate of 
our proposed arrival traffic is determined as follows. Let fc denote the vector defined as: 

J^^^MV^- 



The arrival rate is then given by : 



A = ^^(fi(l — e)f* + t'ie(r* + fc)) = i7+ efc. 

i = 

Lemma 1. There exist queue lengths Qi,Q-2t ■ - Q m, where M is the number of edges in Lmw, such that they 
satisfy: 

QjC^ = QkCki + QiCik = QkCk + Cfc - CkCki > QiCi + c,c,fc. (26) 
Proof: Consider K such that 

= Qk = ^±i±fii^. (27) 
Then, to satisfy equality 2, we need to choose 
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Figure 7: A network with two multiuser links. 



Now, for inequality (3) to hold, substituting for Q, we have, 

^ (CfcCifc + CiCki C/Cfcl C; - „ , 

Qfe^ > [cl- CkCu] + 

Cik Cik 

c,c,k - cl + CkCki- (29) 

From our assumption that Cki/c^ + Ci^/ci > 1, the left hand side of ( |29] l is positive. Hence for ( |29] l to hold, we 
need 

This must be satisfied for every (fc, I) G 



Let p = max 



C-iCki + C^Cik C[Ck 

Then choose ii' such that (iC + CfcC^, — c\)/ck > p for all links (fc, I) e Lmw Then, (|27] | gives the values of Qj 

and Qfc. The value of Q, corresponding to Q^. is then obtained from ( |28] l. Thus one can find a > and hence 

a positive queue vector Qg satisfying relation (ISTT l. In a similar manner one can find positive Cj, c^, and c, that 

satisfy relation dZST l. ■ 

In the example given below, we find the initial queue state vector Qo for a given Lmw To illustrate the proof, 
we provide the following simulation example: 

We consider the network graph described in Example 6. The individual link capacities are 10 packets, while 
the fixed rate for the BC links (1,4) and (2,3) are (5,8) and (8,5) respectively. Considering Lj/h- to be the set 
{(1, 4), (2, 3)}, and by choosing = 0.5 * 20/19 * (0.8 [10 10 0] + 0.6 [0 10 0] + 0.5 [0 10]) and j7 = 0.5 * 
[8 6] + 0.5 * [0 5 8 0], we obtain a"^ = 19/20, implying that the efficiency 7 < 19/20 for this network graph. 
We then simulate the unstable traffic pattern described in the proof of Theorem 2 and plot the queue sizes in Fig. IB] 

The figure shows the queue sizes increasing over the observed time slots, thus verifying Theorem 2 for Example 
6. Packets arrive according to rate allocation vector [8 6] with probability 0.5, and the vector [5 8] with 
probabiUty 0.5. Additional packets may arrive into the queues with probability 0.0002. The additional packets are 
such that they satisfy the constraint in eq.(24), in the proof of Theorem 2. 




Appendix C 
Proof of theorem[3] 

The proof idea can be summarised as follows: Every candidate MW subset Lmw of a tree network graph Q has 
a node da, such that every rate allocation vector on Lmw schedules one of the links connected to do- Any two 
convex combinations of rate vectors will therefore satisfy the multiuser local pooling conditions for one of these 
links. We prove this by using a linear program to represent the multiuser local pooling conditions. 

Lemma 2. For the candidate MW subset Lmw, let a\t'™ denote the highest value of CTm for which Lmw satisfies 
the multiuser local pooling conditions specified in Definition |i] M is defined as a matrix whose columns consist of 
all rate allocation vectors in the set ^lmw- ^ denotes the matrix whose columns are all the rate allocation vectors 
in the set ^^lmw- ^l^o, let e and e denote all ones column vectors of length \^lmw\ '^^'^ \^^mw\ respectively. 
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ctm"'^ can then be represented as the solution to the linear program given below: 

max T (30) 

subject to xM" ^ e 
f M° >z T. 

— O5 if {ki G Lmw and k ^ Lmw^ 
where 



M°. 



M,,.c,, + M,,.c,fe, if{k,l)€L 

MWl 

M,j otherwise. 

otherwise. 

where M^j denotes the index of the matrix M. 

Proof: CTm*^^ can be written as the solution to the following hnear program which is obtained from the 
multiuser local pooling conditions. 

inf Um (31) 

<T,3,/3>0 

subject to (T {Md). ^ (M/5)j, for point-to-point link j, 

CT (Ma), c,, + (Ma), c,, ^ (M/3),c,, + (M/3),c,, 

for multiuser Unks (fc, Z), 

anda'e = l, /3'e = l,a>0, ^ ^ 0. 

Here (Ma), denotes index I of the vector Ma. Setting 7 = era, and denoting 1^ as the edge paired with k in the 
multiuser hnk (fc, Z), the dual function of the hnear program is given by: 

inf [a + y x,[{m)i-{'^l)i] + 

712 

^x4(M/3),c,,, + (M^),,c,,, - (M7),c,,, - (M7),,c,,.] 

l! = l 

+ y(7'e-CT) + ^(/3'e-l)|. 

Since 

(M7),= ^ (M,,)70 and (M7), = (M«)7., 
the dual problem can be expressed as 
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Figure 9: Example of a possible Lj^m- for a tree network. Here L^n- — {I2, (^1.^2), ('3, '4), le, {h, h)}- The multiuser 
links are the pairs of edges enclosed by ellipses in the figure. 



max inf < ct(1 — y) — z + ye'-f — > > a;iMij7j 

^>lO,y,s (yn,B I — ^ — ^ 

j = l i=l 



j = i ^ fc=i ' 
'^^A'/;v' / "2 \ N 

3=1 ^ fc=i ^ ^ 

The above maximization can then be reduced to 

max —z 
subject to 

^ i = l 

f^x,[{M,^c,,^+M,^^c,j\ >0, j-1,...,|r.„^i 

fc=l ^ 

i=i 

■n-2 
fc=l 

The dual problem in (|30] | then follows by defining t = — z, and extending the length of x to |f tjjj,,, | so that it has 
\£L,^.ny I — fii — "-2 zero elements. ■ 
The dual problem is to find an x ^ that maximizes the value of t for which the constraints in (l30l l are satisfied. 
From the dual problem, the set Lmw satisfies 1 -local pooling if r = 1 is a solution to the problem, which implies 
that there exists an x, such that constraint ( l30l ) is satisfied with equality. 

We now show that for every Lmw, one can find an x to satisfy the equality constraint in ( l30t . Since Lmw is a set 
of links from a tree network graph, it satisfies one or both of the following conditions: 

(1) Lmw has an isolated point-to-point link, i.e., it consists of two nodes of degree one each that are only connected 
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to each other. 
(2) Lmw has at least one node with degree 1. 
If L 

Mw has an isolated point-to-point link, let us denote the link by I. We then set the Z*'* index of the vector x, 
corresponding to link I to be one and all other elements of x to be zero. Since link I, being an isolated link is 
always served under the node exclusive interference model, every column vector of M^^^^^. has its Z"* index equal 
to one. Such an x yields r = 1 for the dual problem in ( l30l l. and hence Lmw satisfies 1-local pooling. 
When -La/,, satisfies condition 2, we focus on the node to which the node of degree 1 is connected and denote it 
by do. An example of Lj,,„ and the node d,, is shown in Fig. |9] We now describe the construction of the vector 
X. The elements of this vector depends on the type of links in connected to do. The links can belong to 

following types: 

(1) point-to-point links connected to dg. If j is a point-to-point link connecting d,,, then we set Xj = 1/cj. 

(2) Multiuser links that have node dg as transmitter or receiver If a multiuser link {m,n) G L^iv is such that 
edges m or n are connected to do, and if the point-to-point links m and n are not in Lmw, then we set x^ — 

and a;„ = e^„+c|„/ 

(3) If a multiuser link (m, n) e Lmw is such that edges m or n are connected to dg, and if say point-to-point links 
m e Lmw, then set a;„ = ^1 — ^ j /c^^. Note that for each m, the assignment is possible for only one element 
x„ from our assumption that Y„^ n F„ = {"I>}Vm, n e £. 

Finally we set all other indices in x to zero. We now show that if x is chosen as defined, it satisfies the constraints 
in ( l30l l for t = 1. We first note that under the node-exclusive interference model, amongst the links that share do 
as a common node, only one of the links may be active. As a result, every column j of M" satisfies one of the 
following conditions: 

(i) M"^ — Ci, for only one point-to-point link i G Lmw that is connected to do, in which case M° = for any any 
other link k G Lmw, k ^ j having do as one of its nodes. 

(ii) M°^j — c^n, M"^. — c^^^ + c^^ for a multiuser link {m,n) G Lmw where at least one of edges m or n are 
connected to do. In this case M"^ = for any other edge I G £lmw> ^ '^o as one of its nodes. 

(iii) M^^ = c,„, M°j — c^c„^„ for a multiuser link, such that {m,n) G Lmw and n G Lmw, where m has do as 
one of its nodes. Again, M"^ = for any other edge I G £lmw^ ^ ^'^^ ^'^^ °f ^'■^ nodes. 

Additionally, we observe for M" that if condition (iii) is satisfied by some column j for some link (m, n) G Lmw, 
then any column that allocates non-zero rates to m, n satisfies condition (iii). This is because MGMW always 
selects link m as it gives priority to the point-to-point link over the multiuser link. In a similar manner, we can 
conclude that if for some column j, condition (ii) is satisfied by some link (m, n) G Lmw then every column that 
allocates non-zero rates to link (to, n) satisfies (ii). This is a consequence of the fact that condition (ii) implies that 
(TO,n) G Lmw while m ^ Lmw, n ^ Lmw- 

We now exploit the structure of M° to show that x'M" — e is satisfied by the x that we have constructed. Consider 
any column j in M". If j satisfies condition (i) for some i, then Xi — l/c; yields the inner product of x and column 
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j to be 1. If j satisfies (ii) for some {m,n then a;„, = 0, a;„ = l/(c^„ + c^^J again yielding the inner product 
as one. A similar conclusion holds when j satisfies (iii) for some {m,n) with x„^ = l/c„„ x„ = 0. Since every 
column j satisfies one of the three conditions, x'M° = e is indeed satisfied. It now remains to be verified that 

x'M" < e. 

Consider any multiuser link {m,n) G L^w such that the edges m or n have do as a node. Then for any column 
i, (M°^^, M" J may take on the values (c„, c„c„„), (0, c„c„m), and (c,„ri, cl^^ + c^„). When a;„, — l/c„, a;„ = 0, 
the inner product is less than or equal to one. If — 0, x„ — l/(c^„„ + c^„,J, then the inner product is less than 
or equal to one since cj;^^ + c^^^ > max{c„c„^, c^c^„. If column j has a non zero rate for a point-to-point link 
connected to do, the inner product is equal to one. 

our construction of the vector x ensures that the constraints in ( l30b are satisfied for r = 1. This implies that 1 -local 
pooling is satisfied by L^w Since every Lmw satisfies 1 -local pooling, from Theorem [T] MGMW is throughput 
optimal for network graph Q. □ 
As a corollary of |2l we have the following result that yields a lower bound on the value of crf/^"'. 

. ll^'^ll^ n?^ 

^^-»'-max,,,.,,,_,||r-T||,' ^ 
This result follows by first setting all the non-zero elements of x to 1, i.e, by setting x^ = 1, for every j satisfying 
j G Lmw, or G Lmw for some multiuser link One can then normalize x by max^gi |r^^^^^ , ||r^T||i to 

obtain relation ( [32] l. 

Appendix D 
Proof of Theorem|4] 

Proof: We first show that if the multiuser local pooling factor of a graph (?(V,f) is cr^/' then the network is 
stable under the variable rate MGMW algorithm for all arrival rate vectors A satisfying A G cr*^A. 
We consider the fluid limit model of the system as defined in the proof of Theorem [T] Consider the times when the 
derivative exists for all / G £. Let L„{t) denote the set of links with the highest weight at time t. The weight 

of each link in Lo{t), given by eq. ^ is denoted by w„^{t). We show in Lemma [3] that w„^{t) is differentiable 
almost everywhere. The derivative of the weights of links in Lo{t) is then given by 

^(7j (t)cj TO is a point-to-point link j, 

w^{t) = \ A {q,{t)c:Xt) + qdt)c:At)) rn is a 

multiuser link {k,l). 

Let L{t) denote the set of links from Lg{t), which have the maximum derivative of the weights, 

L{t) = {to G C I w„,(t) = max Wi{t)}. 

ieto(t) 

Then, there exists a small 5 such that in the interval {t,t + S), links in L{t) will have the highest weights in {t, t + 5), 



I.e., 



min w^{t) > max w^(t). (33) 

meL(t) meC\L{t) 
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Let c(r) e Ci(t) be the rate assignments under which the links in L{t) achieve the maximum weight at times 
T G {t,t + S). Since variable rate MGMW picks links from L{t) first, the rate allocation vector selected by variable 
rate MGMW, when projected on L{t) will be an element of R^^lt]- We show in Lemma |4] that the service rate vector 
7r(i), under variable rate MGMW, when projected on L{t) is a convex combination of the elements of R^*!) ■ Using 
Lemma 3, along with the fact that A lies strictly within cr'^A., we show in Lemma |5] that there exists m e L{t) 
such that 

w„,{t) < -e*, (34) 

where e* > 0. Since all links in L{t) have the same derivative, ( [34] i implies that < — e*, Vm G L(t). Now, we 
can consider the following Lyapunov function V{t) :— max^^^ w,„. For V{t) > 0, we have that 

— V(t) < maxi&„ < (35) 

dt+ '"6C 

where j^V{t) ^ lim, 

4,0 is the right hand derivative of V{t). This implies that the largest weight must 

decrease in the time interval {t,t + S). Since the above inequality holds almost everywhere in t, the negative drift 
of the Lyapunov function implies that the fluid limit model of the system is stable and hence by the result in ifTTI . 
the original system is also stable. Thus the efficiency of variable rate MGMW is atleast as large as the multiuser 
local pooling factor. ■ 

Lemma 3. For all m G L{t), the link weight w^{t) is differentiable almost everywhere in t>Q. 

Proof: Suppose m is a point to point link j. The weight w^{t) = qj{t)cj is absolutely continuous since 
Cj is a scalar constant and qj{t) is absolutely continuous. When to is a BC link (fc, /), we use the Implicit 
Function Theorem ifTSl to show that w,„{t) is absolutely continuous. Since w^{t) maximizes the inner product 
{{Cki,Cik) , ((7fc(t), qi{t))), the optimal rate pair can be expressed as ri{t) {qk{t), qi{t)), where r]{t) is the proportion- 
ality factor. Using the Implicit Function Theorem, we can express 77 as a continuously differentiable function of q^. 
and Qi. Since qk{t) and qi{t) are differentiable almost everywhere, and by writing w^{t) as ql{t)rj{t) + qf{t)rj{t), 
we obtain that w„^{t) is also differentiable almost everywhere. ■ 

Lemma 4. The service rate vector Tr{t) projected on L{t) is a convex combination of the elements of R^^^ly 

Proof: We know from Eq. (l33T l that min„gi,(t) w^nir) > ma.x„f,c\L{t) w{t) in {t,t + S). This implies there 
exists an rii such that for all n > Ui, 

min w"{t) > max w"{t). 

ieL(f) " pe\L(f) 

Let r"(T) be the rate allocation vector chosen by variable rate MGMW in time slot t G {t,t + S). Then for 
all n > rii, r^{T) projected on L{t) belongs to the set R^L^t]^ where c(t) G Ci(4) is the transmission rate vector 
associated with the links in L{t) at time r. Consider the service rate vector tt" defined in terms of the cumulative 
service process as: 

S:{n{t + 5)) -S^[nt) ^^^^^^^ 
'no 

We now show that tt" can be expressed as a convex combination of the rate allocation vectors in R^'t). We first 
consider an edge I, such that {k, I) is a multiuser link in L{t). From Eq. ([36]), tt" can be expressed in terms of the 
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rate allocation vector in each time slot in the interval {t,t + S), so that 



„ ^ C'+'-rc-)^- (37) 



Considering only those time slots when r, (r) serves link {k, I), and denoting as the fraction of the time I is 



served as a point-to-point link, we have 

^ 



nS 



oiiCi (38) 
+ aic,. (39) 



+ a,Ci+a.,c;^[t), (40) 



n5 

where — '=° ^' . Since c*^(t) is continuous in t > 0, there exists an e such that |c*^(Ti) — c*j.(t)| < e. Then, 

for all n > rii, we can bound the first term of Eq. (|40l ) as 

E -=0 /„" [c;, (t) - (t)](ir efc max. r_ 

; J — - ^- ^'^^^ 

no no 

Since c*^.{t) is continuous in Vn > rii, as (5 —> 0, we have e 0. From Eq. ( |40l i, we note that by making e 
arbitrarily small, tt" can be expressed in terms of the transmission rates in c{t) alone. Also, since the transmission 
rate of any point-to-point link k is fixed for any c(r), tt^ can be expressed in terms of c{t). Hence, for any edge I 
in L{t), we can write tt" for n > ni as, 

TT" ^ V^-— (42) 

nd 

where r" (r) e 't, . The second term in Eq. (l42l i is a convex combination of the rate allocation vectors in R^'^'j, . 
As tt" — 7- 7r,(f), e — >^ 0, and hence n{t) is a convex combination of the rate vectors in Rl*!)- ■ 

Lemma 5. Given that A G cr'^A, f/iere ejciifi a link m G L(t) ^mc/; that < — e*, /or iome e* > 0. 

Proof: If A G o-l,A, A projected on the subset L{t) is of the form al^fi, where /i is a convex combination 
of the rate allocation vectors in Ri,(t). Then, using Lemma |4] and the fact that L(t) satisfies (7*^ local pooling, we 
know that there exists a link m G L{t) satisfying Eq. (|44| ). given by: 

Xj < TTj (t) m is point-to-point link j, or (43) 

Afc < TTk and A, < TT, (i) m is a multiuser link (fc, Z). (44) 



Eq. ( l44l l implies that there exists a link m G such that 
for point-to-point link j, or 



^^^^c^ = {\it)-7fM)c,<0, (45) 



= (^.(OcLlO + %it)c:,{t)) < 0, (46) 
for multiuser link (k, I). Eq. (l45T l follows if Eq. ( |43] | is satisfied. We then provide the following argument to show 
that Eq. (l44l i implies Eq. ( |46] | for a multiuser link (fc, Z): Suppose ( l46l l is not satisfied and w„ > 0. Since Eq. (l44l i 
implies that ^ < and < 0, there exists a (5,, where S, > 0, such that 
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q,{t + S,) <q,{t), q,{t + 5,) <q,{t), and (47) 
w„,{t + S,)>w„,{t). (48) 
Let (c*,(i + Si),c;^{t + Si)) = argniax^^^^^^^^^^ q^:{t + 61)0^1 + q,{t + 6i)ci^. Then, Eq. gill implies that 

q,{t)c:Xt + (50 + q>{t)cl{t + d,) > wUt + 5i)- 
However, from our definition of weight q^{t)cl^{t+ 5 ^)+ qi{t)c*^{t+ 5 < w„^{t), and therefore w^{t) > w„^{t+5^), 
which contradicts Eq. (l48T l. Hence, ( |46] | must be satisfied. Then, by defining e* > similar manner as done in Eq. 
(fTsT l. we conclude that w„^{t) < — e*. ■ 

Appendix E 
Proof of Theorem[5] 

Suppose, for a candidate MW subset Lmh- with associated transmission rate vector c° e Ci^^^^^., there exists a 
(Jm > 0, and a pair of vectors iJ, fl satisfying (TmIJ. > v,; where is a convex combination of rate vectors in Ri°,,^^ 
, and /2 is a convex combination of rate vectors in ^L^i^^r- Then, in a manner similar to the proof of Theorem |2] we 
can construct a traffic pattern with arrival rate \ = v + ek, where e > is arbitrary and fc ^ is fixed, such that 
the system is unstable under the variable rate MGMW policy. The arrival rate A is arbitrarily close to, but outside 
the boundary of the region ctmT. Let v — X^^'j^^ Q!;?^;, and X^^^^ cti ~ 1, where G ^Vf^ny- We now specify the 
statistics of the arrival pattern in each time slot: Let the initial queue lengths Q(0) = 0. At each time slot, t >1, 
we pick a vector G with probability a^. Conditioned on r^, one of two events may occur. 

1. With probability 1 — e, 

(a) Cj packets arrive into the queue of point-to-point link j, Vj G Lmw such that ri{i) = Cj. 

(b) c°(fc) and c°(Z) packets aiTive into the queues of every multiuser fink {k,l) G Lmw for which ri{k) > 
and r,(0 > 0. 

2. With probability e, 

(a) Cj + Cj packets arrive into the queue of every point-to-point link j G L^w for which rj(j) — Cj. 

(b) Cf^ + c°(fc) and c, + c°{l) packets arrive into the queues of every multiuser link {k,l) G L^w for which 
r,(fc) > and r,(/) > 0. 

(c) Cj packets arrive into the queues of all edges j G Sl^w ^^^^ that ri{j) = 0. 

Cj, Ct,i, and Ci^ are fixed positive quantities defined for every point-to-point link j G Lmw, and (fc, /) G Lmw- They 
satisfy the following conditions: 

CjC, = CkC''{k) + CiC''{l), \fj, {k, I), G Lmw, 
Q C°(fc) = C°(0 C„ V(fc, G Lmw. (49) 

In the following lemma we show that these quantities exist. 

Lemma 6. There exist fixed quantities Cj > 0, corresponding to every edge j G £lmw' ^^^^ ^'^'^^ ^'^^J satisfy the 
conditions in ( |491 l. 
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Proof: The proof of Lemma |6] is similar to the proof of Lemma 1 and follows by equating the quantities in 
( |49] l to a constant K and then computing each Cj, c^i and Cj^ in terms of K. ■ 
We now discuss the behaviour of the queues in Lmw in Lemmas 6 and 7. As a consequence of Lemmas 6 and 
7, we will show that the queues become unstable under the described arrival traffic. 

Lemma 7. Suppose the queues in L^,^y satisfy the condition in ( |49l l as stated below: 

QjCj = QkC''{k) + Q,c°{l) V j, {k, I) e L^„, 

Q,c°{k)=c"{l)Q„ \^{k,l)eLMw (50) 

Then, if the queue length Qj, for every j € ''^ increased to Qj + Cj, the new queue lengths will again satisfy 

relation ( I50I) . 

Proof: From relations ( |49l ) and dSOl l. we observe that (Q^ + Cj)cj ~ {Q^ + Ck)c°{k) + (Q, + c,)c°(Z). Also, 
using dSOb we obtain (Q, + c,) c°(/c) = c°(Z) (Q^ + c^). ■ 

Lemma 8. Suppose at the end of time slot t — 1, the queues in Lj^-^y satisfy the relation <I50I) . Then, at the end of 
time slot t, with probability 1 — e, the length of all queues in L^w remains unchanged, and with probability e, the 
lengths of every queue in Lmw increases by a fixed positive quantity; i.e., if Qj = Qj{t + 1) — Aj{t) denotes the 
queue of j G ^Lf^i^y '^^ ^"'^ of time slot t, then with probability 1 — e, Qj{t) = Qj(t — 1) Vj G ^tj/w" ^'^'^ ^'^^ 
probability e, Q,{t) ^ Q,{t) + c,, £ 

Proof: We consider the arrival statistics of packets that we defined previously. Conditioned on vector being 
picked, packets only arrive into those links that have non zero rates in f*. We now show that the weight of links 
that have a non-zero rate in Ti, when calculated after the arrival of packets, strictly dominates the weight of all 
other links in Lmw', and hence, MGMW picks the rate allocation vector . Consider first the packet arrivals as 
described in event 1. Queues in Lmw, by assumption satisfy relation ( l50b at the end of slot t. At the beginning of 
time slot < + 1, after packets have arrived, the following is true: For all multiuser links [k, I) G Lmw satisfying 
r^{k) > and > 0: 

(c°,c°)= argmax [(Qfc + c°)cfc, + ((5, + c°)c,j,], and (51) 

(Q, + cl)cl + [Q, + c°)c° >Q™C + Q-c:- (52) 

y{m,n) G Lmw satisfying ri{m) = r'i(n) — 0. Since and Qi satisfy ( l50b . (Q^ + c°) and Q, + c°) also satisfy 
dSOt , implying that the queue length ratios remain unchanged. ( fSTT ) follows from the fact that the queue length ratios 
remain unchanged and the pair (c°, c°) lies on the boundary of TZki- ( l52b is obtained by noting that the weight of 
links that are not served in Ti, remains unchanged. In a similar manner, for point-to-point links j G Lmw satisfying 

r.{i) > 0, 

Q,{t) + c,)c, > max (q, (i)c,, + Q„(t)c:, ) (53) 
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for all links I, {m,n) e L^w satisfying ri{l) = ri(m) = ri{n) = 0. Relations ( |52] l and ( |53] ) together imply that 
MGMW schedules rate allocation vector over Lmw in time slot t. Since the number of packets arriving into 
queue j equals r^Q), at the end of time slot t, Qj{t) = Qj{t — I), Vj e ^Lj^t^- 

Next, consider the case when packets arrive as described in event 2. Applying Lemma |2l the queue lengths + 
Cj, Vj e ^ijjH ^^^1 satisfy ( fSOb . Since the additional packets that arrive into each j G equals 7'i(j), the 

same argument used for packet arrivals in event 1 implies that MGMW again picks rate allocation vector r^, and 
hence, Qj{t) = Qj{t — 1) + c,, at the end of time slot t. ■ 
We can now describe the queue evolution pattern for our packet arrival traffic. The initial queue state = 
0,Vj' G Sl,,.,^, satisfies dSOb . From Lemma [8] the queue lengths do not decrease, and the event that each queue 
length increases by a fixed positive quantity occurs infinitely often. Hence, the system is unstable under our proposed 
arrival traffic. The arrival rate of the traffic pattern is then evaluated as follows: A = '^Zllo ^ ^)^» + '^Zlll 
where k{j) ^ c,, Vj e Slj,,^- □ 



